Drag on particles in a nematic suspension by a moving nematic-isotropic interface 
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We report the first clear demonstration of drag on colloidal particles by a moving nematic-isotropic 
interface. The balance of forces explains our observation of periodic, strip-like structures that are 
produced by the movement of these particles. 
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Colloidal dispersions of small particles in nematic liq- 
uid crystals are a novel, interesting type of soft matter. 
The difference from ordinary colloids arises from the ori- 
cntational ordering of the liquid crystal molecules and 
the resulting structure in the colloid. Topological defects 
Jj], [| || and additional long-range forces between the 
colloidal particles Q are immediate consequences of this 
ordering. The nematic-induced interparticle interaction 
brings a new range of effects to the system: supermolec- 
ular structures || 0, B], cellular structures 0, pd| , 
and even a soft solid jl2j can be observed. Colloidal 
dispersions in liquid crystals also have a wide variety of 
potential applications |0j |. 

A range of problems similar to those of polymer dis- 
persed liquid crystals also arise in nematic colloidal dis- 
persions. The nematic ordering makes it difficult to sus- 
pend small particles in a liquid crystal host. Particles of- 
ten segregate into agglomerates distributed nonuniformly 
in the cell. The resulting spatial distribution of the parti- 
cles is difficult to control. Our research explores the fac- 
tors that affect the spatial distribution of these particles 
and indicates ways to control the complex morphology of 
these systems. 

In this paper we report the first demonstration of drag 
on colloidal particles by a moving nematic-isotropic (NI) 
interface. We calculate a critical radius above which the 
particles cannot be captured by the moving interface. 
We predict that this critical radius is sensitive to the 
viscous properties of the host liquid crystal, the value 
of the anchoring coefficient of the liquid crystal on the 
particle surface, and the velocity of the moving interface. 
Most important, we can move particles of specified radius 
and can control the spatial distribution of these particles 
in the cell. 

In order to understand how the particles are moved 
by the nematic-isotropic transition front we used parti- 
cles of different size as well as particles made of different 



materials. In the first part of our experiments we used 
nearly monodisperse spheres of silica (R = 0.005, 0.5, and 
1/j.m). To prove our predictions based on initial experi- 
mental observations and to demonstrate the controllabil- 
ity of the spatial distribution of particles in anisitropic 
colloidal suspensions, we used large polymer particles, 
R = 8/im |l7j]. In all cases particles were dispersed at 
concentrations of qh = 1 — 5wt %, in the liquid crystal 
5CB at room temperature (25°C) (the isotropic-nematic 
transition temperature of pure 5CB is Tni ~ 35°C). The 
sample was subjected to ultrasound in order to uniformly 
disperse the particles in 5CB. Some of the preparations 
were made at higher temperatures, in the nematic or 
isotropic state of the liquid crystal. 

The homogeneous mixture was deposited between two 
polymer-covered glass substrates and heated above the 
nematic-isotropic transition point. The cell thickness 
was 100/im. The homogeneous suspension was observed 
under crossed polarizers using an optical microscope, 
Fig. 0a. 

The mixture was cooled to a temperature below the 
transition point. Depending on the rate of cooling, we 
observed different structures. Fast quenching to room 
temperature (cooling rate 10°C min -1 ) resulted in phase 
separation and formation of a cellular structure, with 
particle-free nematic domains separated by particle-rich 
regions, Fig. |l|b. Properties of these structures have been 
reported previously [ll| |l2). Decreasing the cool- 
ing rate, we observed formation of a striped structure 
(Fig. |l|c). The particle-rich regions were no longer form- 
ing a cellular structure but were arranged in a set of 
stripes, separated by particle-free regions. Using opti- 
cal microscope images we postulate that we have large 
nematic and isotropic domains separated by a moving 
interface. The direction of the stripes is parallel to the 
moving interface (the interface was moving from the left 
to the right of the cell in the geometry depicted in Fig. 111). 
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FIG. 1: Polarized microscopy images of different structures, 
depending on the cooling rate: (a) colloidal particles are 
dispersed homogeneously in the isotropic phase; (b) cellular 
structure: cooling rate I0°C min -1 ; (c) stripes: cooling rate 
0.1°C min , velocity of the interface v w 3fim sec , (d) 
root structure: cooling rate 0.0I°C min -1 , v ~ 0.5(im sec -1 . 
Silica particles, R = 0.005/im, concentration <f> = Iwt %. The 
NI interface is moving from the left to the right. The long 
side of the images is 1mm. 



The spatial period of the striped structure depended on 
the cooling rate, as well as on the particle size. Increas- 
ing the particle size, as well as decreasing the cooling 
rate, resulted in an increase of the spatial period. We 
also noticed that the stripes do not appear if we have 
considerably larger silica particles, R > 0.5/xm. Also, de- 
creasing the cooling rate resulted in a chaotic merging of 
stripes and formation of a "roof -like pattern, Fig. [yd. 
These results indicate that the particles are pushed by 
the moving nematic-isotropic phase transition front. 

Using optical microscopy we were able to directly 
observe the movement of larger (R ~ 8/im) but less 
dense polymer spheres. We could see these particles be- 
ing moved by an advancing nematic to isotropic phase 
boundary. Fig. shows pictures of this moving front 
taken at different times. Clearly the particles are pushed 
by this advancing front, remaining in the isotropic phase. 

A simple analysis accounts for these experimental ob- 
servations. We consider several mechanisms which con- 
tribute to the total drag force acting on a particle at 




FIG. 2: Snapshots of the moving nematic-isotropic interface 
made between parallel polarizers taken with the time interval 
At = 6min. The nematic phase is on the bottom part of 
the picture. Polymer particles, R = 8/j.m. One can observe 
thickening of the particle-rich phase because of the capturing 
of the particles by the interface. Velocity of the interface is 
v ~ 3/Ltm sec -1 . The long side of the images is 1mm. 



a NI interface. First, the surface tension coefficient 
might differ at the particle-nematic or particle-isotropic 
part of the interface. An additional pressure caused by 
the curvature can be given by P = 2a /R, where R is 
the radius of the curvature (in our case it is the ra- 
dius of the particle) and a is the surface tension coef- 
ficient. This pressure contributes to the total drag force 
as F a — 27r(crjv — — (d/R) 2 ], with the amplitude 

growing linearly with the droplet radius. Here d is the 
distance from the particle center to the interface. The 
value of Ao" = un — 07 depends on the surface treat- 
ment of the particles and is unknown (and difficult to 
measure) for our system. However, the order of mag- 
nitude can be estimated from the change of the sur- 
face tension coefficient of the glass-5CB interface [Q, 
Act « 10 -2 - 10 -3 dyne cm -2 . 

Second, the particle creates long-range distortions of 
the director in a nematic phase. To minimize elastic dis- 
tortion energy, the nematic tries to expel particles into 
the isotropic phase. The elastic forces have two origins: 
due to the director deformations in the bulk nematic and 
due to the anchoring of the director at the particle sur- 
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face. An estimation of these contributions can be done 
by dimensional analysis. For the surface contribution, 
the only combination which has dimension of force is 
WR, where W is the anchoring coefficient. Therefore, 
the surface contribution to the drag force is proportional 
to WR, F s — WRg s (d/R), where g s (x) is a dimensionless 
function of the penetration depth, d/R. 

One can have two different situations for the bulk 
contribution. For weak anchoring, WR/K << 1, the 
bulk contribution is proportional to the squared char- 
acteristic deviation of the director, /3q ~ WR/K [ p5| . 
Now W 2 R 2 /K has the dimension of force, yielding Fb = 
W 2 R 2 /Kgb(d/R). In contrast, in the case of strong an- 
choring, WR/K >> 1, the anchoring does not enter the 
elastic contribution, and Ft, — Kg\,{d/R). Here, again, 
gb(x) is a dimensionless function. For 5CB Q, and 
typical values of the anchoring energy, W ~ 1CP 3 — 
10 _4 dyn cm -1 , WR/K << 1 for silica particles and 
WR/K rj 1 for polymer particles. Therefore, for silica 
particles, we have the weak anchoring regime. In con- 
trast, polymer particles provide strong anchoring of the 
director. We also note that, when particles agglomer- 
ate, the effective radius increases and we have a strong 
anchoring regime even for small particles. 

Finally, there is a friction drag contribution, which, in 
the first approximation, is given by the Stokes formula, 



F n = —6irRr]v |16|. The total drag on the particle is the 



sum of all contributions, F, 



drag 



F„ 



Solution of Newton's equations of motion with -Fdrag as 
a force completes the description of the particle dynam- 
ics. It is clear, however, that small heavy particles cannot 
be moved by the interface. The maximal radius can be 
estimated from the conservation of linear momentum. To 
capture a particle of mass to, the interface has to transfer 
to it a linear momentum tow. If we assume that the par- 
ticle does not move (or it moves much slower than the 
interface, which is valid for massive particles) then the 
total linear momentum transferred to the particle reads 



Fdr&gdt 



1 



Fd iag (x)dx, (1) 



Here we assumed that the interface touches the particle 
at time t\ and leaves it at time ti, x — vt. Substituting 
-Fdrag we obtain 



Rn 
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S h W 2 /K 



(2) 



where Si are geometrical constants, p is the density of the 
particle, Ar is the final displacement of the particle due 
to the drag force. 

Several important conclusions can be drawn. First, if 
the particle is too big, the moving interface is not able to 
transfer sufficient linear momentum to it. Only particles 
with R < i? max (v,W,a) will be captured by the inter- 
face. From eqn (0) one can see that i? max v~ 2 , i.e. 
only a slowly moving interface is able to capture the par- 
ticles. The estimate of this velocity is given by the zero of 



the denominator of eqn (||), v W/^/Kp ~ 1mm sec -1 . 
This is of the order of the limiting velocity for the cel- 
lular structure we observed in our experiments: if the 
interface moves more slowly, then stripes appear, other- 
wise the cellular structure forms (see Fig. [I]) . 

The main conclusion is that i? ma x is a function of the 
material parameters, i.e. can be effectively controlled, for 
example, by changing the surface treatment of the par- 
ticles (anchoring energy W). Increase in the anchoring 
energy leads to an increase of i? max - Moreover, strong 
enough anchoring favors formation of a defect near the 
particle |[], ||, contributing to an even higher energetic 
barrier created by elastic forces. 

On the other hand, if the particle is captured by the 
interface, the elastic force scales as R 2 , and the oppos- 
ing viscous drag scales as R. Therefore, there is a min- 
imal radius, i? m i n , starting from which particles will be 
dragged by the interface. If the particle is dragged by the 
interface at a constant speed, then Fdrag = 0, yielding 



Rm'in — 



6in]v — 2nAa — j s W 
lh W 2 /K ' 



(3) 



where -fi = gi(0) are some constants. Eqn (||) implies 
that, in order to be moved by the interface, the particles 
have to be big enough. Only in this case can elastic forces 
overcome viscous drag. Substituting values typical for 
5CB and using the slowest cooling rate, we obtain i? m i n ~ 
O.Ol^m which qualitatively agrees with the minimum size 
of silica particles we were able to move. 

To explain formation of the striped structure, we note 
that, in practice, particles aggregate into clusters. While 
an aggregate moves, it captures more and more particles, 
growing in size. The anchoring parameter WR/K also 
increases and we switch from the weak anchoring to the 
strong anchoring regime. The bulk elastic contribution is 
then proportional to the elastic constant K and the elas- 
tic force is no longer growing as R 2 . Therefore, at some 
R c , the friction drag overcomes the elastic contribution 
and the aggregate breaks through the interface. A stripe 
forms and the particles start to accumulate again. The 
condition Fd rag = gives the critical size of the aggregate 



R c = 



IbK 



6Trrjv — 2nAa 



(4) 



which is about lpm for typical experimental values. 

From conservation of mass one can show that the ra- 
dius of the aggregate increases linearly with time until it 
reaches R r , 



R = R + -<p—vt. 

4 p a 



(5) 



Here i?o is the initial radius of the aggregate, p a is the 
density of the aggregate. Therefore, the distance between 
two stripes is given by 



A w 40" 



■1 Pa n 
PLC 



(6) 
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and is of the order of 0.1mm, again in qualitative agree- 
ment with experiment, Fig. [l|. 

In conclusion, we have demonstrated moving colloidal 
particles by a moving nematic-isotropic interface. We 
have also determined the factors such as particle size, 
anchoring energy and speed of the moving front, that 
control the particle movement. By controlling these fac- 
tors we can control the morphology of the colloid and its 
physical properties. Such control is necessary to develop 



and optimize these colloids for specific applications. 
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